r 

AD-A152  926  A  MATHEMATICAL  MODEL  FOR  CALCULATING  DETECTION 

PROBABILITY*  OF  A  DIFFUSION  TARGET(U)  NAVAL  POSTGRADUATI 
SCHOOL  HONTEREV  CA  M  SISLIOGLU  SEP  84 

UNCLASSIFIED  F/G  12/1 

:■  1 

1 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

1 

L 

■ 

k. 

m 

A 152  026 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


$ 


DT1C 

ELECTE 
APR  4  1985 


THESIS 


Approved  for  public  release;  distribution  unlimited. 


85  03  18  115 


SECURITY  CLASSIFICATION  OF  '''his  pace  (When  Date  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

(mm 

4.  Title  (end  Subtitle) 

A  Mathematical  Model  for  Calculating 
Detection  Probability  of  a  Diffusion 
Target 

5.  TYPE  OF  REPORT  4  PERIOD  COVEREO 

Master's  Thesis; 
September  1984 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR^ 

Mucahit  Sislioglu 

8.  CONTRACT  OR  GRANT  NUMBER/' »J 

S.  PERFORMING  ORGANIZATION  name  ANO  AOORESS 

Naval  Postgraduate  School 

Monterey,  California  93943 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 

11.  CONTROLLING  OFF|CE  NAME  AND  ADDRESS 

Naval  Postgraduate  School 

Monterey,  California  93943 

12.  REPORT  DATE 

September  1984 

13  number  of  pages 

69 

U  MONITORING  AGENCY  NAME  4  ADORES S(l(  different  from  Controlling  Office) 

15.  SECURITY  C l.  ASS.  (of  this  report) 

is*  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 

_ J 

IS.  DISTRIBUTION  STATEMENT  (of  rM  s  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  .'of  the  .  -.fr«cf  entered  In  Block  20,  if  different  from  Report ) 

18.  SUPPLEMENTARY  NOTES 

19  KEY  WORDS  ( Continue  on  reverse  side  if  neceesmry  end  Identity  by  6/ocA  nuirb*r) 

Diffusion  target,  random  tour  model,  detection  probability 

20  ABSTRACT  f Continue  on  reverse  side  If  neceseery  end  Identify  by  block  number ) 

The  primary  objective  of  this  stretchy-  is  to  derive  a  mathe¬ 
matical  model  to  predict  the  detection  probability  of  a  target 
which  moves  randomly,  according  to  a  two-dimensional  diffusion 
model.  This  model  assumes  that  there  is  a  stationary  searcher 
which  has  a  " cook ie-cutter ,r  sensor  with  radius  R.  In  order 
to  construct  this  model,  a  Monte  Carlo  simulation  program  is 
used  to  generate  detection  probabilities.  It  is  demonstrated 

DD  1  JAN  73  1473  COITION  OF  1  NOV  «S  IS  OBSOLETE 


MO;*  LF-  01  4-  660! 


SECURITY  CLASSIFICATION  OF  this  R*GE  Omte  Entered) 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  Dmtm  Ertffd) 


Block  20  Contd. 

that  this  model  can  be  used  asymptotically  to  predict  an  upper 
bound  detection  probability  of  an  '“equivalent"'  random  tour 
target.  .■>  ,  ,  7 '  •  '  Z-i  / 1 '  [  . 


S  N  0102-  LF-  014-  6601 


2 


SECURITY  CLASSIFICATION  OF  this  PAGE (Whtt  Dtm  Enf*r*rf) 


Approved  for  public  release;  distrinution  unlimited. 


Author: 


Approved 


A  Mathematical  Model  for  Calculating 
Detection  Probability 
of  a  Diffusion  Target 


by 


Mucahit  Sislioglu 

Lieutenant  Junior  Grade,  Turkish  Navy 
B.S.,  Turkish  Naval  Academy,  1978 


Submitted  in  partial  fulfillment  of  th» 
requirements  ror  the  degree  of 

MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 

from  the 


NAVAL  POSTGRADUATE  SCHOOL 
September  19  84 


vCv-uAt.  h 

- Kneire-T7”flarsfi'aTI7  - 

Dean  of  Information  and  Policy  5 ciences 


ABSTRACT 


The  primary  objective  of  this  study  is  to  derive  a  math¬ 
ematical  model  to  predict  the  detection  probaaility  of  a 
target  which  moves  randomly,  according  to  a  two-dimensional 
diffusion  model.  This  model  assumes  that  there  is  a 
stationary  searcher  which  has  a  "cookie-cutter"  sensor  with 
radius  R.  In  order  to  construct  this  model,  a  Monte  Carlo 
simulation  program  is  used  to  generate  detection  probabili¬ 
ties.  It  is  demonstrated  that  this  model  can  be  used 
asymptotically  to  predict  an  upper  bound  detection  prob¬ 
ability  of  an  "equivalent"  random  tour  target. 
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I.  DESCRIPTION  OF  THE  DIFFUSION  MODEL 


A.  INTRODUCTION 

The  main  objective  of  this  thesis  is  to  finl  and  test  an 
experimental  mathematical  model  whicn  predicts  the  proa- 
anility  of  detecting  a  two-dimensional  tacget  by  a 
stationary  searcher.  This  model  will  be  shown  to  provide  an 
upper  bound  for  the  probability  of  detection  by  a  stationary 
searcher  of  a  target  conducting  a  "random  tour"  'Ref.  1] 

B.  DESCRIPTION  OF  DIFFUSION  MODEL 

1.  The  Searcher  location 

The  searcher  is  assumed  to  be  located  in  the  center 
of  a  sguare  search  region  of  area  A.  This  location  is  held 
fixed  during  the  search  period.  The  searcher  has  a  detec¬ 
tion  capability  over  a  disk  of  radius  E.  The  detection 
probability  of  a  target  inside  of  this  disk  is  1  and  outside 
is  0.  The  searcher  thus  has  a  "cookie—  cutter"  sensor  with 
detection  range  F..  [fief.  2  ] 

2.  The  Target  Starting  Position 

The  target’s  starting  position  is  uniformly 
distributed  over  the  sguare  search  region  A. 

3.  Motion  of  the  Target 

In  our  diffusion  modti,  the  target  moves  randomly 
over  the  area  A  as  a  diffusing  particle  which  ceflects  off 
ti.e  area  boundaries.  The  diffusion  constant  is  D ,  which  has 
dimensions  of  area  per  unit  of  time.  In  any  time  interval 
of  length  At  that  does  not  contain  a  boundary  reflection. 


components  of  the  target's  position  on  the  X  and  Y  axes 
suffer  increments  which  are  independent  of  all  previous 
increments  and  which  are  eacn  distributed  normally  with  mean 
0  and  variance  D£t. 

Still  ignoring  boundary  effects,  the  diffusion 
assumption  results  in  the  target's  location  at  time  t  having 
a  circular  bivariate  normal  probability  distribution  with 
mean  of  the  starting  position  and  variance  of  Dt. 

Thus  the  probability  density  of  the  tacget's  loca¬ 
tion  at  time  t  is 


£  (x,y,t)  = 


27T  Dt 


exp  - 


fx-u,)2  » 

2  Dt 


(i.i) 


where  (Ux#Uy)  is  the  target's  starting  position.  Adding 
the  effects  of  boundary  reflection  significantly  complicates 
the  calculation  of  f(x,y,t)  and  leads  to  the  necessity  of 
using  simulation  to  attack  this  problem, 

4.  Detection 

Detection  occurs  whenever  the  target  enters  the 
searcher  circular  detection  disk  which  has  a  radius  R. 


C.  DIFFUSION  SI  MU  L ATIQN  MODEL  (DIFSIM) 


A  Monte  Carlo  simulation  computer  model  (DIFSIM)  is  used 
to  genarate  detection  probabilities  for  this  diffusion 
model.  This  program  is  written  in  FORTRAN  and  designed  for 
use  at  the  Naval  Postgraduate  School ( NPGS) .  It  uses  the  new 
version  of  the  NPSS  Random  Number  Genarator  Package,  called 
LLRANDOMII  in  order  to  generate  Uniform  and  Normal  random 
numbers. 
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1  *  Inputs 

•  Area  size.  A,  in  square  nautical  miles. 

•  Diffusion  constant,  D,  in  square  nautical  miles  per 
hour. 

•  Padius  of  detection  disk,  R,  in  nautical  miles. 

•  Number  of  replications. 

•  Detection  period  as  an  hour. 

•  lime  increment,  ^t,  for  each  discrete  step  in  minutes. 
2.  function inq  ox  Program 

The  initial  target  position  is  s^le^ted  from  a 
bivariate  uniform  distribution  over  the  search  region  A. 
Subsequent  target  positions  are  determined  by  a  discrete 
approximation  of  the  diffusion.  We  make  the  following  defi¬ 
nitions, 

X=x  component  of  current  location 
Y=y  component  of  current  location 
X,  =  x  component  of  new  location  at  the  end 
of  time  increment  At 

Y’  =  v  component  of  new  location  at  the  end 
of  time  increment  At 

Then , 

X’=X  +  €^(D£t)'/4 
Y  ,  =  Y+  G^DAt)'^ 

where  ^  and  6^  are.  drawn  independently  from  a  standard  normal 
distribution. 

In  this  model  a  5  minutest  is  used.  Different  time 
increments,  varying  from  1  minute  up  to  15  minutes,  have 


been  tested  and  5  minutes  has  been  accepted  as  a  good  value. 
For  smaller  time  increments  the  simulation  program  took  too 
long  in  computer  execution  time.  As  shown  in  Figure  1.1 
there  is  no  significant  difference  in  probaoiiity  curves 
between  5  minute  increments  and  smaller  time  increments. 
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When  the  target  er.cour.ters  a  boundary,  a  reflection 
is  uade  to  Keep  the  target  inside  tre  search  area.  In  a 
target's  Y’  position  after  a  reflection  is  giver,  as  follows: 

Y' <  0  =>  Y*  becomes  -Y ' 

Y’>  a  =>  Y'  becomes  2a-Y ' 

where  a  is  the  length  of  a  side  of  tne  square  search  area  A. 
The  target  reflects  ir.  the  X  direction  in  a  similar  manner. 

Detection  occurs  whenever  the  target  enters  the 
detection  disk.  This  event  can  be  defined  analytically  as 
follows : 

(X--J-)  2  f  (Y~~)  2<  R2 
3 .  Cut  put 

For  each  time  t,  the  simulation  output  is  the  ratio 
Nq/Nt  where 

N0=number  of  replications  giving  a  detection  by  time  t, 
UT=total  number  of  iionte  Carlo  Replications  used  in  the 
simulation. 
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Figure  1.2  Diffusion  Model 


II.  CONNECTION  BETWEEN  R ANDOH  TOOR  AND  DIFFUSION  MODEL 


A.  DESCRIPTION  OF  RANDOM  TOOR  MODEL 

In  the  random  tour  model  considered  here,  the  target  is 
assumed  to  move  at  a  constant  speed  and  to  make  course 
changes  at  random  times.  Each  new  course  is  drawn  2  cor 
uniform  distribution  on  [0,2JT]-  The  lengths  of  the  time 
between  course  cnanges  is  exponentially  distributed  with 
mean  1/A- 

An  analytic  expression  for  the  probability  density  of 
the  target’s  position  after  a  random  tour  of  length  t  was 
derived  in  [Ref.  1]  Given  the  target's  initial  position  at 
the  origin  of  a  two-dimensional  coordinate  system,  this 
expression  is 


g(x,y,t)=-2m--jl 


™”-exp  t  ( 1-/1 -5*;)  +exp''*'t  |S(S" 1  j 


(2.  1) 


V=target  speed  (nautical  mile  in  per  hour) 
A=course  change  rate  (hrs  ) 
t=time  (hrs) 


The  Dirac  £_function  component  of  g(x,y,t)  arises  from 
the  fact  that  with  probability  e  the  target  does  not 
finish  the  first  step  by  time  t.  In  other  words,  the  target 
makes  no  course  change  through  time  t.  Therefore,  its  prob¬ 
ability  mass  is  concentrated  on  the  boundary  of  a  disk  of 
radius  vt  and  centered  at  the  origin. 
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,  ^ii'orr  tie  ii.t  c  j  i :  o:.  i  r.  (2.2),  *  u 

V  t  i> 


A  t 


e  x  n  f  A  t  J 1-x  ■  (  d  x 


1  “  x 4 


0 

=  X  t  [  (1-U2)  exp  ^  Atuj 


d  u 


Xt  , 
=  e  - 1 


Xt 

1  <V  ^  d. 


(Xt) 


2  e 


A  t 


(Xt-l  +  e  ^  C )  -  1  (2.4) 


(Xt) 


From  (2.2)  and  (2.4)  we  then  have 


i 

(2.5) 

» 


E[R2]  =  2  — (Xt  -  1  +  e  Xt) 

X2 

2 V 2  t  /1_  1 -e~ 

X  \  Xt 


In  the  calculation  of  E£ E2  ]  for  diffusion  model  with 
diffusion  constant  D  ,  the  target's  initial  fosition  is 
assumed  to  he  at  the  origin  of  a  two-dimensional  coordinate  | 

system.  V.e  have 

E[H2]=E[X2+Y2  ]=E[X2  ]+E[Y2  ]  (2.6) 


I 


Since  X  and  Y  are  maepead  ent  ana  uncorrela ted,  they  have 
normal  marginal  distributions.  So, 

X— ii  (0 , Dt ) 

Y~N  (0  ,Dt) 


E[  X2]=Var[  X  >  ( 3£  X  ])  2=Dt+0=Dt 
Z[  Y2  ]=  Var  [  Y  ]♦  (S[Y])  2=Dt+0  =  Dt 

If  we  substitute  these  E[X2]  and  EJ  Y2  ]  values  into  (2.6),  we 


E[  E2  ]  =  e[  X2  ]+E[  Y  2  ]=Dt  +  D  t=2Dt 


(2.7) 


As  described  in  [Ref.  3]  for  tae  random  tour  model,  as  t 
goes  to  infinity  the  Central  Limit  Theorem  requires  that 
9 (X/Y/t)  becomes  asymptotically  circular  bivariate  normal 
with  mean  p  =0  and  variance  o"2=7 t/fi.  This  result  can  be 

obtained  by  using  the  formula  (2.5)  and  letting  t  go  to 
inf i nity . 


2V* h  I  -  e~ht  2V*t 

Lim  E[R2]=Lim  - ^ - (1 - ^ - )= - ^ — 

t - >  OO  t - >  <=>© 


(2.6) 


By  comparing  the  equations  (2.7)  and  (2.8),  it  is  seen 
that  as  t  becomes  large,  a  random  tour  can  be  approximated 
by  a  diffusion  with  diffusion  constant 


(2.9) 


To  examine  the  relationship  between  a  random  tour  and 
its  ’'equivalent"  diffusion,  two  simulations  were  used. 
Example  results  of  these  two  simulation  programs,  DIFSIM  and 


PASS  are  displayed  m  Figure  2.1.  3oth  programs  are  Monte 
Carlo  search  simulations.  For  PASS  (Passive  Acoustic 
Submarine  Simulation),  [Ref.  4]  the  target  motion  model  is  a 
random  tour.  In  DIFSIM  (Diffusion  Simulation),  the  target 
moves  as  a  diffusing  particle.  In  both  cases,  tne  searcher 
is  stationary  at  the  center  of  a  100  r.m  x  100  nm  search 
area  and  has  a  cookie_c.it  ter  sensor  with  detection  range 
15nm.  The  target  reflects  orf  the  area  boundaries.  Ihe 
initial  target  position  was  selected  uniformly  over  the 
search  area,  and  each  replication  of  the  simulation  was 
allowed  to  continue  until  the  target  moved  within  distance 
15  nm  of  the  searcher. 

To  generate  the  results  shown  in  Figure  2.  1  ,  the 
following  5  different  pairs  of  h  and  V  are  used  as  a  rate  of 
course  change  and  speed  of  target  in  the  PASS  simulation 
model . 


15.  0 


0.25/,  \  1.0/  ,  \  2.25J  ,  \  4.0 


20.0 


25.0 
i  o.25, 


nm/hr 

1/nr 


If  we  use  equation  (2.9),  we  may  get  an  "eguival  ent "  diffu¬ 
sion  constant  103  no2/hr,  by  using  eac.h  dirferent  (A, 7) 
pair.  Thus  100  nm2/hr  is  used  as  a  diffusion  constant  in 
tne  DIFSIM  in  order  to  get  an  "equivalent"  diffusion  model. 

As  demonstrated,  detection  probabilities  ara  asymptoti¬ 
cally  very  close  to  each  other  as  t  increases.  But  during 
the  early  search  hours,  detection  pro baoilities  f or  a  diffu¬ 
sion  model  are  higher  than  the  pr ooabiiities  which  are 
generated  by  the  random  tour  model. 

If  we  recall  the  the  equations  (2.5)  and  (2.3)  we  will 
see  that 


Thus  the  approximation 
model  by  using  the  diffu 
OVER ES1I HATE  of  E[  E  2  ]. 

In  the  diffusion  model 
will  move  a  greater  distan 
target  in  the  "equivalent" 
fore  reasonable  to  expect 
encounter  a  stationary  sea 
target  conducting  an  "equiv 
ture  has  been  supported  by 
supported  is  the  fact  that 
caiiy  identical. 

An  experimental  analyti 
the  diffusion  model  in  the 


of  E[  R  2  ]  for  the  random  tour 
sion  model  always  leads  to  an 

we  may  expect  that  the  target 
ce  from  tne  origin  taan  does  the 
random  tour  modex.  It  is  there- 
that  the  diffusing  target  will 
rcher  more  quickly  than  will  .a 
alent"  random  tour.  This  conjec- 
further  simulation  tasting;  also 
the  two  processes  are  asymptoti¬ 


cal  model  will  be  constructed  for 
next  section. 


Ill,  ESTABLISHMENT  of  the  analytical  model 


In  this  cl  ap  tc  r  ,  and  with  si  muia  ti  on  results  from  TIFSIM, 
an  experiment al  analytical  model  will  be  constructed  for 
predicting  the  probability  of  a  diffusion  target  entering  a 
stationary  disk  by  seme  time  t. 

A,  ASSUMPTIONS 

The  following  assumptions  are  used  in  our  model 

•The  searcher  is  fixed  at  the  center  of  a  sguare  search 
region  of  area  A, 

•The  searcher  has  a  detection  capability  ovec  a  disK  of 
radius  E  with  a  detection  probability  of  1  within  the 
disk  and  0  outside.  (i.e,  a  cookie-cutter  detector). 

•The  target’s  starting  position  is  uniformly  distributed 
over  a  search  region  A. 

•The  target  moves  randomly  over  the  area  A  as  a  diffusion 
particle  with  diffusion  constant  0, 

•  The  target  reflects  off  the  area  boundaries. 

•A  target  can  be  detected  only  once  by  the  seaccher 

B.  CLASSIFICATION  OF  THE  VARIABLES 
The  variables  in  our  model  are, 

•First  detection  probability,  P. 

•  Search  area,  A  (nm2)  . 

•Searcher  detection  disk  radius,  R  (am). 

•  Target  diffusion  constant,  D  (nm2/hr)  . 
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•  Time ,  t  (hr)  . 

The  first  detection  probability,  ?,  is  the  dependent  vari¬ 
able.  The  remaining  variables  A ,  h ,  D  and  t  are  independent. 
That  is. 

?=  f(A,H,D,t) 

C.  CCHSIEUCTION  DF  THE  MODEL 

Figure  3.1  shows  four  p^ots  of  the  proDaoility  of  a 
target  detection  by  time  t  as  estimated  by  the  Monte  Carlo 
simulation  DIFSIM.  If  we  look  at  tnese  carves,  we  will 
observe  that  all  of  them  have  an  increasing  trend  and  tday 
approach  1  asymptotically.  It  also  appears  as  if  the  second 
derivative  must  be  negative  everywhere.  Figure  3.2  plots  1 
minus  the  same  data  cn  a  log  scale.  The  fact  that  these 
plots  are  very  nearly  linear  suggests  the  following  func¬ 
tional  form  P  (t)  . 

P(t)  =  1-cxe"^  (3.  1) 

where  a.  and  p  are  determined  by  tne  problem  parameters  R,  A 
and  D.  After  conducting  23  separate  simulations  with 
different  values  of  R,A  and  D,  the  author  is  convinced  that 
the  form  of  P(t)  is  approximately  exponential.  This  thesis 
attempts  to  establish  values  of  <x  and  as  function  of  E,A 

and  D,  to  allow  equation  (3.1)  to  be  a  reasonable  estimate. 


l  gure 


Figure  3.2  {-  P  (t)  Plottou  on  Logarithnic  Scale 
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•  Submodel  for  ot 

It.  our  model,  we  know  that  the  target  starting  posi¬ 
tion  is  uniformly  distributed  over  the  search  area  A,  and 
the  searcher  has  a  detection  capability  over  a  disk  w.i ic.1 
covers  an  area  Hp2  with  probabilty  1.  So,  ws  may  expect 
that  at  the  beyinning  of  the  search,  i.e.,  when  t=0,  detec¬ 
tion  probability  will  be  egual  to  1TE2/A. 

If  we  subsititute  t  =0  in  equation  (3.1),  we  get 


P(0)=  1-ex 

Then 

TT  R1 

E(t=0)  = — — - =  1-  cx 

which  implies  that 


(3.2) 


(3.3) 


2.  Submodel  for  p 

This  submodel  will  include  ail  independent  variables 
E,D,A  and  t,  and  will  therefore  be  more  complex .  We  will 
study  each  independent  variable  separately  in  order  to 
simplify  the  problem  (We  will  change  one  variable  while- 
holding  the  others  fixed.)  The  relationship  between  p  an  1 
these  variables  will  be  estimated.  Eventually  we  will 

combine  these  submodels  for  a  final  submodel. 

a.  The  Relationship  Between  Diffusing  Constant  D 
and  p 

For  this  case,  area  size  A  and  ralius  P  w=r-e 
held  fixed  at  1  0000  nm2  10  nm.  respectively.  Diffusion 
constant  D  was  varied  between  20  nm2/ar  and  303  nn2/hr  The 
simulation  results  are  displayed  in  Figures  3.3. 


Figure  3.3  P(t)  for  5  Diffusion  Constants 
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Lffusion  Constant 


h.  The  Relationship  Between  Area  Size  A  and  |3 

This  time  the  diffusion  constant  D  aid  detection 
radius  R  were  held  at  100  nm2/hr  and  10  nm  respectively. 
The  area  size  A  was  varied  between  lOOOnm2  and  20000  r.m2. 
The  simulation  results  and  the  least  square  estimation 
results  for  the  log  transformed  data  values  are  d  isplayed  in 
Figure  3.5  and  and  Table  II 

Figure  3.6  shows  a  plot  of  the  best  fit  p 
against  area  size  A.  These  points  fit  very  closely  with  the 
power  function 


rea 


TABLE  II 

Exponential  curve  Fitting  For  Dif 
Area  Sizes  R=10  nm  D=100  nn2 


c.  The  Relationship  3etween  Detection  Radius  R  and 

P 

Ajain,  the  same  procedures  were  applied  for  this 
case.  R  was  varied  between  1  and  30  nm,  while  D  and  A  were 
held  fixed  at  100  nm2/hr  andlOOOO  na2  respectively. 

The  simulation  results  and  least  square  estima¬ 
tion  results  are  displayed  in  Figure  3.7  and  Table  III.  In 
Figure  3.8,  the  scatter  plot  of  estimated  J3  values  shows  a 
linear  relationship  between  R  and  p. 

Least  square  estimation  for  this  lin2  is 

ft  =-0.00105  +  0.  00278R  (3.  6* 

pH  0 . 00278 R 


Figure  3*7  P(t)  for  5  Radius  Value 
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Figure  3.8  Best  Fit  Q  vs.  Radius 


We  can  summarize  these  ooservations  is  follows: 


P  o-  D 
O*.  A  “*•* 

P~  H 

where  means  "is  proportional  to"  wnich  suggests  that 


P=K 


RD 

A,s 


(3.7) 


for  the  proper  value  of  K.  To  estimate  K  is  the  final  model 
building  step. 

d.  Estimation  of  the  Coefficient  K 


Ws  can  calculate  the  K  value  for  each  simulation 
with  the  expression 


K= 


P  A 


I'S 


RD 


(3.3) 


where  ji  is  the  "best  fit"  value  for  that  simulation  run. 
Then  with  these  sample  '{,  values  we  may  find  best  overall 
estimate. 


In 

25  additional 
of  31  sample 
table  values 


addition  to  56  simulations 
simulations  were  conducted 
K  values.  The  histogram 
for  this  data  are  displayed 


already  completed, 
to  produce  a  total 
and  the  statistical 
in  Figure  3.9. 
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Figure  3.9  Observed  Values  of 


If  vt  recall  the  ijuatioa  (3,1),  <  j  .  3 )  .a!  (3.7) 

l’U)  =-  1  -  .t-' 


ex 


=  1 


7lR‘ 


A 


A 

Substituting  these  <x  and  £  values  in  equation 
(3.1),  we  derived  our  final  analytical  model  for  first 
detection  probabilities  as  follows 

r  2  RD  r 

P(t)  =  1  -  (1  -  - Zf—)  e  K  A  i  .5  0.9) 

where  ?  (t)  refers-  to  the  probability  of  first  detection 
occurs  on  or  before  tine  t. 

D.  VERIFICATIOH  OF  THE  MODEL 

1  .  Dimension  Analysis 

From  equation  (3.9)  we  see  that 


K  R  D  t 
A's 


must  be  dimensionless.  This  implies  that  the  coefficient  K 
must  be  dimensionless.  (If  we  had  set  the  power  of  A  to 
1.49  versus  1.5,  then  the  dimension  of  K  would  be  nm0*2,  not 
a  natural  unit.) 
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2.  Sensitivity  Analysis  for  Independent  Variables 

If  ve  hold  fixed  all  independent  variables  E,D,t  and 
change  the  area  size  A,  we  will  observe  that  as  we  increase 
the  A,  the  probability  of  detection  decreases  and  vice 
versa.  This  result  is  demonstrated  in  Figure  3.  10.  As  we 
increase  the  size  of  the  search  area,  more  area  will  be 
available  for  the  target  to  escape  frotr.  the  searcher. 
Therefore  we  may  expect  lower  detection  probabilities  for 
larger  search  area  sizes. 

A  similar  sensitivity  analysis  is  applied  for  tne 
other  independent  variables  R,D  and  t.  If  we  look  at  the 
results  displayed  in  Figure  3.10,  we  may  observe  that  as  we 
increase  these  variables,  the  detection  probabilities 
increase  simultaneously.  These  results  seem  reasonable, 
because  as  we  increase  the  searcn  time  or  detect  ion  radius, 
we  may  have  more  chanses  to  detect  the  target.  Also  an 
increase  in  the  D  value  means  that  the  target  will  travel 
more  distance  during  any  time  interval  and  will  thus  be  more 
likely  to  enter  the  detection  disk. 


Figure  3.10  Sensitivity  Analysis  for 


3.  Final  Verification 

There  was  no  actual  data  available  fron  real  life 
observations.  Therefore,  the  output  of  DIFSIM  is  used  for 
final  verification  of  the  model.  For  this  purpose,  combina¬ 
tions  of  the  following  independent  variables  were  used  both 
in  our  analytical  model  and  as  input  to  DIFSIK. 

D=  40,  80,  140,  200  Nm*2/Hr 

B=  2.  5,  1  0,  20,  25  Nm 

A=  4000,  8000,  12000,  16000  Nm*2 

The  outputs  are  displayed  in  Figure  3.  .11,  3.12, 

3.13.  It  is  observed  that  the  simulation  and  model  prob¬ 
abilities  are  generally  very  close  tc  each  other.  Only 
during  early  search  hours  do  the  simulation  curves  sometimes 
go  above  the  model  curve.  This  means  that  our  model 
predicts  fewer  detection  than  the  diffusion  simulation  model 
during  the  early  search  hours.  Since  this  difference  is  at 
maximum  .03  or  .04,  we  conclude  that  the  model  provides  a 
good  fit,  which  gets  better  for  larger  time  t.  For  small 
t,the  model  appears  to  underestimate  the  proa  ability  of 
detection. 
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DIFF.  CONS.-  40  Nm»2Alr  £  f  OIFF.  CONS.-  00  Nm*2/Hr 


Model  Verification  Results 


Results  (A  Varies 


01 FF.  CONS.-  100  Nm*2/Ur  ~  I  OtFF.  CONS.-  100  Nm* 2/Mr 


Pesults  (R  Varies 


7  m  t i  i  <  -  :  I  t  <  d  ,  t  i  i.:  *  -  . :  *  r  :  t  < 

t  ■  •  •  ■  j  - 1 .  i  ;  i ;  i  c  '  '  r.  T 


Z:  »-  fact  tiiut  the  II?  fir  rcobebiiity  oi  ctt».c::Dn  i  y  tine  t 
g-.  rm:  rally  e;;cjec"  that  predicted  ry  t:i«  model  that 

the  ir.stan  tan  eo  us  probability  of  detection  produced  by 
DlrSI^  exceeds  (at  least  in  early  hours  of  the  simula¬ 
tion).  To  test  this  hypothesis,  DIFSIfl  was  run  viti:  the 
parameters 

D=  1  00  nm2/hr 

E= 1 0  Dd 

A  =  1  0000  nm2 


The:,  for  each  5  hour  period  between  0  and  00  hours,  the 
least  squares  Lest  fit  was  obtained.  (That  is,  the  best 
fitting 


1  -  (1 


-  3 1 
e 


v:as  calculated).  These  values  are  plotted  in  Figure  3.14, 
ar.d  appear  to  approach  from  above  the  model  value  of 

24, 7  x  100  x  10  =  q247 

10000 

Thus,  in  this  case  at  least,  it  appears  that  the  instanta¬ 
neous  xrobafility  of  detection  starts  at  some  high  value  and 
decreases  asymptotically  to  a  steady  state  value  given  by 
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tii a  model  presented  here.  This  observation  remains  tc  be 
proved  in  general. 
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ns tdiitdjieous  Dotection  Fate  (j3)  vs.  Ti:ae 


iv.  rr  CD/.riLisiic  a n  *. i y :* i s  of 


^odf: 


CUluIAIIVE  DllLCHOU  PL  01  ABILITY  rUlCIIGN 

Since  the  following  ;>ro:K  rtits  for  general  cu;.jLiUv- 
.  r  ;» i  l:  Lili  t  y  function  hold  for  our  :  -ticr:l  rod*_l,  ^u-i- 

tiuii  (3.5),  we  nu  y  assure  tut  t.hx£  oiel  also  nr  iu  -ji.u  a 
emulative  density  function  (cdf)  for  first  detection  t-^.e 

t. 

These  properties  are 
1.  Lin  F (t)  =  1 


t  —  >  °° 
Lin  F  (t)  = 

t  —  >  CO 


-K^ 
e  A  r 


1 


2.  F  (t)  is 

If  we 
r espect 


a  nondecr easing  function 

take  the  the  first  derivate  of 
to  t,  we  get 


F(t>  with 


ar  ( t) 

dt 


0 


since  this  equation  always  has  norregutiva  values,  we 
nay  assure  that  F (t)  is  a  nondecr casing  function. 

3. r  ( t )  is  a  continuous  function. 


Therefore  we  nay  define  cdf  as  following. 


T. 


Since  tnis  integration  is  not  egual  to  1/  the  equation  (4.2) 
doesn’t  represent  a  proper  density  function  for  t. 

C.  EXPECTED  FIRST  DETECTION  TIME 

Let  F (t)  denote  the  cumulative  nondetection  probability 
function  (cndf )  . 

F  (t)  =  P  (No  detection  up  to  time  t)  =  p  (T>t) 


=  1-P(T«t)=1-?<t) 


The  following  formula  can  be  used  to  find  the  expected 
detection  time  E£T]. 

E[  T  ]=  Jf(t)  dt  (4.4) 

o 

If  we  substitute  (4.3)  in  (4.4) 

emr|r(|-iS£)e-K7S<]d+ 

n  *■  •* 


t 
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f  01 


Exfectec  first  detection  tines  are  displayed 
different  diffusion  constants  p )  ,  area  s^zes  (Ai  and  detec¬ 
tion  radiuses  ( F: >  in  Figure  4.  1. 

D.  CONDITIONAL  DETECTION  PROBABILITY  FUNCTIONS 

If  we  assume  that  there  will  he  no  detection  at  the 
beginning  of  the  search  period,  we  may  derive  tae  following 
condi tional  cdf . 

F0  (t)  =E  (Det  ec  t  ion  up  to  time  t/no  det.  at  time  0) 

zP[T4i/T>0];  PCT>0,nU 

prr>oj 

P[o<T£i] 

P[T>Oj 


If  we  substitute  t=0  in  equation  (4.3),  we  jet 
P  (T  >0 )  =  F  (0  )  =  ( 1 — LaC_)  (4'7) 

and 


P(0<T$t)  =  F  (t)  - P  (t=0) 


By  usiny  (4.7) 
equation  (4.6) 


F„(0=  |  _ 


a  nd 
,  we 


(4.8) 

have 
ft  O 
A,.S 


as 


a  dominator  and  numerator  xr.  the 


(4  3) 
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This  function  (4.9)  rapreseuts  a  cumulative  probability 
function  of  ar.  exponential  distribution  witti  pararete: 

RO 

K — ^7-j - .  By  using  this  fact,  conditional  expected  first 

detection  time  E[T0]  can  be  defined  as  follows: 


,lT 


(4.  10) 


E[T0]=E[r/No  detection  at  tine  0  ]=  — 

Also,  we  car.  write  conditional  detection  probability  der.sit; 
function  f0  (t)  in  the  following  form. 


e‘"  a- 


n) 


If  we  compare  equations  (4.5)  and  (4.10)  ,  we  will 
observe  that, 

AyfTft  /  A,s 

KRO  KO  v  KRO 


so  E[T]<E[T0].  This  inequality  means  that  tha  conditional 
first  detection  time  is  greater  than  the  unconditional  first 
detection  time.  We  can  gat  this  conclusion  intuitively  by 
thinking  that  we  have  an  opportunity  to  detect  the  target  at 
the  beginning  of  the  search  period  in  the  unconditional 
case.  Therefore,  for  the  unconditional  case,  w?  may  ex;ect 
an  earlier  first  detection  time  than  would  be  possible  for 
tha  conditional  case. 
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(oOOll/iTlN)  IW1SN03  ftOlSOJ 


» 


Figure  4.2 


Conparison  of  F  (t)  and  F0(t). 
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7.  THE  RESULTS  AND  APPLICATIONS  OF  THE  MODEL 

A.  EXPECTED  DETECTION  TIME  FOR  RANDOM  TOUR  MODEL 

As  we  showed  in  Chapter  II,  a  diffusion  tar  jet  jives  an 
upper  fccund  detection  probability  for  the  "equivalent" 
random  tour  model.  Therefore,  the  diffusion  model  expected 
detection  time,  which  is  estimated  by  equation  (4.5)  should 
be  a  lower  bound  for  random  tour  expected  detection  time. 
That  is,  for 

E[ T0]=Dif fusion  Model  Expected  Detection  Time,  and 

E[ Tg  ]=  Fandom  Tour  Model  Expected  Detection  Time,  we 
have 


B.  ONE -DIMENSIONAL  DIFFUSION  MODEL 

In  this  thesis,  two-dimensional  diffusion  motion  was  the 
basis  for  our  model.  The  exponential  type  curves  were  used 
to  estimate  this  model’s  outputs.  In  addition  to  this 
model,  the  one- di mensional  diffusion  model  is  simulated  by 
computer  program  DIFSIM1.  In  this  model,  the  target  moves  on 
a  line  segment  L,  according  to  diffusion  constant  D.  The 
target's  starting  position  is  selected  uniformly  over  this 
line  segment  L.  Detection  occurs  whenever  the  target  hits 
the  designated  end  point.  The  target  reflects  off  the  ether 
end  point  of  line  segment. 

The  results  of  different  simulation  results  are 
displayed  in  Figure  5.1.  Exponential  curves,  which  were 
obtained  by  using  the  least  squares  estimation  method,  were 
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used  to  estimate  the  outputs  of  one-dimensiona 1  diffusion 
model  as  we  did  for  tvo-dim ensionai  case. 

So  we  nay  expect  that,  for  three  or  more  dimensional 
diffusion  models,  we  may  use  exponential  type  curves  which 
are  generated  by  a  different  set  of  parameters.  For  three- 
dimensional  case,  these  parameters  can  be  defined  as 
follows : 

V=Volume  of  the  cubical  search  space. 

R=Radius  of  the  cookie-cutter  detection  sphere. 

D=Diffusion  constant. 

T=Detection  time. 

C.  APPLICATIONS  OF  THE  MODEL 

Our  model  can  also  be  used  to  estimate  the  final  detec¬ 
tion  probability  of  a  system  which  includes  more  than  one 
independent  sensor.  As  an  example,  we  may  use. the  following 
scenario : 

Ke  want  to  use  n  sonobouys  in  order  to  detect  a 
diffusing  target  in  an  area  A.  Each  sonobouy  has  a 
cookie_cutter  detection  capability  over  a  disk  with  radius  K; 
.  Each  sonobouy  will  be  independetly  located  on  the  center 
of  a  square  subsearch  area  A; and  operated  for  a  time  period 
t.  If  we  make  the  following  assumptions,  we  may  use  equa¬ 
tion  (3.9)  to  estimate  the  overall  detection  probability  of 
this  sonobouy  pattern  at  the  end  of  the  search  period  t. 

-  (£*•-) 1/1 

13  | 

n 

i*i 

where  R  is  the  effective  detection  radius. 


Figure  5.1  One-Dimensional  Diffusion  Model 


APPENDIX  a 

DIFSIH  C3MP0TER  PROGRAM 

In  order  to  give  access  to  tae  logic  used  in  buildin^ 
the  simulation  models  DIFSIM  and  PIFSI.11,  a  compLete  program 
listing  is  included  in  this  Appendix  following  the  list  of 
variables  used  in  the  simulation  models. 

LIST  OF  VARIABLES 

RE?  =Number  of  replications. 

MAX  =Detection  period  as  a  minute. 

E  =Radius  of  detection  disk  in  nautical  miles. 

DIF  =Diffusion  constant  D  in  square  nautical  miles  per 
hour. 

SIDE=The  length  of  the  scuare  search  area  side  in 
nautical  miles. 

AREA=Area  size  A  in  sguare  nautical  miles. 

INC  =Time  increment  t  for  each  discrete  step  in 
minute  s. 

PROB=Probability  of  detection 
POSX=X  component  of  target  position. 

PGSY=Y  component  of  target  position. 

DI5T=The  distance  between  the  target  location  and  the 
center  of  the  detection  disk. 

1X1  =Seed  number  for  uniform  random  number. 

1x2  =Seed  number  for  standard  normal  random  variable. 
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